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Abstract: The hadron spectrum of SU(2) lattice gauge theory with two flavours of 
Wilson quark is studied on an 8^ x 16 lattice using all-to-all propagators, with particular 
emphasis on the dependence on quark chemical potential fi. As fi is increased from zero 
the diquark states with non-zero baryon number B respond as expected, while states 
with B = remain unaffected until the onset of non-zero baryon density at /x = 771,^/2. 
Post onset the pi-meson mass increases in accordance with chiral perturbation theory 
while the rho becomes lighter. In the diquark sector a Goldstone state associated with 
a superfluid ground state can be identified. A further consequence of superfiuidity is 
an approximate degeneracy between mesons and baryons with the same spacetime and 
isospin quantum numbers. Finally we find tentative evidence for the binding of states 
with kaon quantum numbers within the baryonic medium. 
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1 Introduction 

At large baryon chemical potential fiB the properties of QCD are expected to change as 
the system moves from a confined nuclear matter phase to a deconfined quark matter 
phase where the relevant degrees of freedom are quarks and gluons. At low temperature 
Tand high /i^, the attraction between quarks is believed to be sufficient to promote di- 
quark Cooper pairing leading to a colour superconducting ground state. Weak-coupling 
techniques can be used at asymptotic densities and have revealed a superconducting 
phase known as the colour-fiavour locked phase. However as density is reduced towards 
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phenomenologically reasonable values, the precise nature of the ground state appears 
very sensitive both to additional parameters such as isospin chemical potential and 
strange quark mass, and also to the nature of the non-perturbative assumptions made 
in the calculation. It seems natural to use Lattice QCD to investigate these issues, but 
unfortunately whilst the lattice has been used very successfully to investigate QCD with 
T > 0, the well known "Sign Problem" has made progress for fis/T ^ 1 impossible. 

Orthodox simulation techniques can be applied, however, to the case of two colour 
QCD (QC2D) with gauge group SU(2). Whilst this theory differs in important ways 
from QCD, for instance in having bosonic baryons in the spectrum, and in having a 
superfluid, rather than superconducting, ground state at large /i_B, it remains the simplest 
gauge theory in which a systematic non-perturbative treatment of a baryonic medium 
is possible. Moreover, recent simulations [1] have provided evidence that there exist two 
distinct forms of two color matter: the dilute Bose gas formed from diquark bound states 
which forms at onset, ie. for /i^ > /iso = M^^, in which superfluidity arises via Bose- 
Einstein condensation of scalar diquarks, is supplanted at larger densities by a deconfined 
"quark matter" phase in which a system of degenerate quarks is disrupted by BCS 
condensation at the Fermi surface. Studies in this regime may have qualitative or even 
quantitative relevance for QCD quark matter, particularly in the non-Goldstone sector. 
For instance, Schafer [2] has stressed how the impact of instantons on the excitation 
spectrum at high baryon density could be elucidated by lattice simulations. 

In this Letter we study the /i^-dependence of the hadron spectrum in both meson 
and baryon sectors of QC2D with Nf = 2 flavours of Wilson quark, which we find to 
vary dramatically as the onset from vacuum to a ground state with non-zero baryon 
density is traversed. We build on the pioneering work of the Hiroshima group [3j, who 
found that beyond onset the pion mass did not change noticeably, but the rho meson 
became significantly lighter, so that the level-ordering is reversed. Two baryon channels 
were also studied but no significant /i^-dependence found. This work was performed 
on small lattices, and only looked at a few states. One of our aims is to improve and 
update their results. We also study the nature of the Goldstone mode associated with 
superfluidity, as done for QC2D with staggered lattice fermions in and expose the 
specifically two colour phenomenon of "meson-baryon" mixing in the superfluid state, 
whereby since B is no longer conserved, states interpolated by mesonic operators qq and 
diquark operators qq have identical quantum numbers and hence exhibit an approximate 
degeneracy. Finally, we make the first measurements of strange meson masses in a 
baryonic medium, and find tentative evidence for bound states of kaons in nuclei. 

2 Formulation 

The gauge-invariant lattice action with Nf = 2 degenerate fermion flavours is [1] 

S = AM{fi)A + ^2M(/x)V^2 - '^j{^iK^2 - ^2K^i), (1) 
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with M the conventional Wilson fermion matrix (with lattice spacing a = 1) 

MM = ^-y-^^J2 - 7.)e'^'^°t/,(x)5^,,+p + (1 + 7.)e-^'-°f/J(y)5^,.-p] , (2) 

K the hopping parameter, ^ the quark chemical potential, and j the strength of an 
SU(2)i®SU(2)/j-invariant diquark source term needed to regularise IR fluctuations in 
the superfluid phase, which should be extrapolated to zero to reach the physical limit. 
The subscript on the fermion fields is a flavour index. The anti-unitary operator K = 
= (^75X2, where C-^^C^^ = — 7J = —7* and the Pauli matrix T2 acts on colour 
indices. A useful relation is 

M^(/i) = -K-f,M{-fi)K-f5. (3) 

The hadronic states examined in this paper are qq mesons and qq, qq diquark baryons 
and anti-baryons. In all cases we use local interpolating operators of the form 'ilj{x)rilj{x), 
ip'^{x)KTilj{x), ip{x)KV%l)'^ [x). The matrix V = 7or'''7o determines the spacetime quan- 
tum numbers of the hadron, with inclusion of the K factor ensuring that mesons and 
baryons with the same F have the same J^. In this letter we will focus on states with 
r G {1, 75, 7j, «757j} with j = 1, . . . , 3 corresponding to G {0+, 0~, 1~, 1"*"}. 



2.1 Fermion Propagators 

The fermion action ([1]) can be written in the form '^A4{fi, j)'$ where ^ = {ipi, ipj, %l)2-, ipj) 
and \l/ = {■ipi,'ijj2 ,ip2,'4'i)'^ ■ It has the form: 
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and A{j) = A{~j). Now consider the propagator 
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where the zero entries arise from the assumption of isospin symmetry. This symmetry 
also implies that 
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where the subscripts denote "normal" and "anomalous" propagation. Anomalous prop- 
agation arises from particle-hole mixing in a superfluid ground state; on a finite volume 
Sa vanishes in the limit j — » 0. 



3 



2.2 Mesons 



The isovector (/ = 1) meson operators are given by il)iTip2-, il'2^'4'i and {ipiTipi — 
%l)2Til)2) / \/2. The charged meson correlator is then 

= - Ti[SNiy,x)TSNix,y)T] + Tr[S'A(z/, a;)rS'^(x, ?/)r^] 

For the neutral meson correlator, the "disconnected" parts made up from the product 
of two traces cancel because of isospin symmetry. The connected parts are: 

- {iJ2ix)TMx)My)^My)), + {Mx)TMx)My)^My)), (8) 

= - Tr[5'iv(?/, x)TSn{x, y)T] + Tr[5'^(?/, x)TSAix, y)T^] 

and so the degeneracy between neutral and charged isovector mesons is manifest. The 
isoscalar (/ = 0) meson = {ipiTipi + '?/'2rV'2)/V2 has the correlator 

(MO(x)MOt(y)) = 2TT[SN{x,x)T]TT[S^{y,y)T] 

- T^'^[SNiy,x)TSNix,y)T] - TT[SAiy, x)TSAix,y)T'^]. 

In this case there is a disconnected term, and the anomalous term has the opposite sign 
to that of the isovector meson. From now on the isovector 0~ and 1^ mesons will be 
referred to as the pion and rho respectively. Note that for j = /i = the difference 
between / = and / = 1 is entirely due to disconnected diagrams, which must therefore 
account for the n - rj' mass splitting. For fj,, j ^ this need not be the case. 

In this Letter we for the first time use lattice techniques to study the spectrum of 
kaons in a dense baryonic medium. We model the kaon states by isovector mesons ([7]) in 
which one of the quark flavours does not feel the chemical potential (for details of how 
this is achieved see Sec. E]). The resulting propagator thus approximates the behaviour 
of an s quark in a medium made up from u and d quarks, although at this exploratory 
stage we make no attempt to give the strange quarks a realistic mass but simply assume 
Ks = The states studied here correspond to K"^ {J^ = 0~) and K*"^ {J^ = 



2.3 Diquarks 

An isoscalar diquark is given by 

D\x) = {ijl{x)KTMx) - ^Ux)KTMx)) (10) 
implying a correlator: 

2{D'^ix)D%y)) = {i,,{x)TK^^{x)i;l{y)KTi;2iy)) - {Mx)rKij^{x)i;^iy)KTi;,{y)) 
- {Mx)rKijl{x)ijUy)KrMy)) + {Mx)TK^p^{x)iJUy)KTiJ^{y)) 
= 4 Tr[5A(x, x)TK] Tr [5a(z/, y)KT] 

+ 2 TT[S^iy, x)TKSNix, y)f ^iT] + 2 Tr[^^(y, x)TKSNix, y)KT]. 

(11) 
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The last two terms arising from connected diagrams cancel if KT^K = T, which is the 
case for T = i'y^'jj, and add up if KT^K = — f implying F G {l,75,7j}. Thus the only 
baryon states available in the 1 = sector are O"*", 0~ and 1~; this is just the Pauli 
exclusion principle in action. The isovector diquark correlator is 

{D'^{x)D\y)) = TT[Sr,{y, x)rKSN{x, y)T^K] - Tr[5^(y, x)rKSN{x, y)KT]. (12) 

Here the disconnected diagrams cancel, and the connected terms only add up for F = 
i757j implying = 1+. 

Spontaneous symmetry breaking of the U{1)b symmetry ip i— > e*"'?/', ip i— >■ -i/^e"*" in 
the superfluid phase (the symmetry is explicitly broken when j ^ 0) can be probed 
by studying two special diquark states in the isoscalar 0"*" sector which we refer to as 
"Higgs" and "Goldstone" [5]. The operators, which by construction yield correlation 
functions symmetric under Euclidean time reversal, are given by 

D±(x) = ^^{x)K^l{x) ± ^l{x)K^2{x). (13) 

where the Higgs is '+' and the Goldstone '— '. After Wick contraction the correlators 
are 

{D\x)D^{y)) = TT[KSA{x,x)]TT[KSA{y,y)] + TT[KSA{x,x)]TT[KSA{y,y)] 

± {TT[KSA{x,x)]TT[KSA{y,y)r]) + TT[KSA{x,x)]TT[KSA{y,y)]) 
- {Ti[SA{y,x)KSA{x,y)K]+TT[SA{y,x)KSA{x,y)K]) 
± {Tr[SNiy,x)KSNix,y)K]+Ti[SNiy,x)KSNix,y)K]) . 

With the choice of diquark source specified in ([T]), in the superfluid phase the U(l)s 
symmetry is preferentially broken by a condensate (D^) ^ 0, and as the name implies, 
in the j — > limit a massless Goldstone mode is interpolated by D~ . The degeneracy 
between Higgs and Goldstone is principally broken by connected diagrams formed from 
normal propagators. 



3 Numerical Method 

We studied an ensemble of gauge configurations generated on a 8^ x 16 lattice at various 
values of /i using a Hybrid Monte Carlo algorithm with the fermion action ([T]) supple- 
mented by a standard Wilson gauge action [1]. The parameters were /3 = 1.7, k = 0.178; 
studies of the string tension suggest a "physical" lattice spacing a = 0.26(l)fm|l and 
studies of the fj, = meson spectrum yield M^a = 0.79(1) and M^^/Mp = 0.80(1) [6]. 
For the most part the diquark source ja = 0.04, though for a few values of fi the series 
ja = 0.06, 0.04, 0.02 was studied in order to permit a j — >• extrapolation. The values 
of /i studied are given in Table [H Configurations are separated by 4 HMC trajectories 
of typical mean length 0.5, and are not in general independent. 

■^This value is based on string tension measurements on a 12"^ x 24 lattice and supplants those 
reported in [TJ [5] . 
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/ia 


ja 


No. of Configurations 


jia 


ja 


No. of Configurations 


0.00 


0.04 


210 


0.50 


0.02 


133 


0.10 


0.04 


200 


0.50 


0.04 


248 


0.25 


0.04 


117 


0.50 


0.06 


76 


0.30 


0.02 


78 


0.55 


0.04 


410 


0.30 


0.04 


133 


0.60 


0.04 


441 


0.30 


0.06 


70 


0.70 


0.02 


178 


0.35 


0.04 


124 


0.70 


0.04 


382 


0.40 


0.04 


175 


0.70 


0.06 


70 


0.45 


0.04 


281 


0.90 


0.04 


65 



Table 1: Number of configuations examined at each value of chemical potential. 



Quark propagators on each configuration were calculated by all-to-all techniques [7] 
using time, spin and flavour dilution. Whilst this equates to constructing 16 x 4 x 2 = 128 
dilution vectors per propagator and thus performing 256 inversions per configuration, at 
/X 7^ this was required for acceptable statistical precision. These were saved to disk and 
stitched into correlation functions according to the formula ([7]), ([9]). ffTTl) . f|T2|) and f|T^ . 
Disconnected diagrams relevant for each state were calculated and saved seperately; with 
the current level of statistics these contributions are both noisy and compatable with 
zero and so only the connected pieces are presented in this initial study. 

To construct kaon correlators, partial quenching was used for the strange propagators, 
i.e. a further set of propagators was calculated on a /z 7^ gauge configuration using 
= 0; j = 0.04). This choice of j ensures degeneracy of eg. K and pion states 
at /X = 0. The two different sets of propagators were then stitched together with the 
relevant operators to form K and K* states. We reiterate that the s quark mass used 
was identical to that of the u/d, and clearly whilst the kaon states contain contributions 
from virtual u/d loops, there is no s-quark sea contribution at this stage. 

To extract masses the meson correlation functions were fitted to a cosh function. The 
time range was adjusted to achieve a stable fit while minimising the obtained x^/^df- 
For yU 7^ states with baryon number B ^ such as the diquarks and the kaons are 
no longer degenerate with their anti-particles. This results in correlators which are no 
longer time-symmetric and must be fitted by a sum of two independent exponentials. 
This can lead to problems with fitting as a 4 parameter fit is more susceptible to noise 
than one with just 2. Naively in the vacuum phase below onset the states' masses 
receive an additive contribution ±^BNc, and so as fi increases the correlation function 
is increasingly dominated by the lighter of the particle-antiparticle pair, which is the one 
with B < 0. For fi > the heavier state rapidly becomes very difficult to fit. 
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4 Results 



All results presented here come purely from analysis of hadron correlators formed from 
connected quark propagators. The result for the pion mass at /x = is M^^a = 0.800(3). 
Chiral perturbation theory (xPT) [8], which models the weakly- interacting dilute Bose 
gas forming at onset, can be used to predict the onset chemical potential /Xo where the 
transition from vacuum to the superfluid phase occurs. This prediction, which strictly 
speaking only applies when there is a separation of scales between the pion and heavier 
hadrons, is that fio = M.„/2. In our case, this suggests that the onset in the j ^ limit 
occurs at fio^ — 0.4. As fi is increased we found that more iterations are required to 
perform an inversion of the quark matrix [Ij, especially once > fio due to the large 
density of small eigenvalues of in the neighbourhood of the origin in this regime. 

In principle all our results should be extrapolated to the "physical" limit j — 0, but 
unfortunately available resources preclude a systematic study for all /x. Here we follow 
[Ij by studying ja = 0.02, . . . , 0.06 at three representative /i points: just below onset, 
just above onset, and well into the superfluid phase. The results for and Mp are 
shown in Fig. [TJ Within the limits of statistical precision a linear extrapolation j — is 
valid. It is notable that the pi and rho masses have opposite slopes as functions of j, and 
that the slopes change sign across the onset transition between /xa = 0.3 and /xa = 0.5. 



|Xfl=0.3 



|j,fl=0.5 



|Xfl=0.7 




0.02 0.04 0.06 0.02 0.04 0.06 0.02 0.04 0.06 



ya 



Figure 1 : Pi and rho masses as a function of diquark source strength j for three representative values 
of /X, together with their values extrapolated to j — 0. 
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Results for the meson spectrum as a function of /i are displayed in Fig. [2l At = 



1.5 



0.5 




_L 



_L 







0.1 



0.2 



0.3 0.4 
[ia 



0.5 



0.6 0.7 



Figure 2: Meson masses as a function of for ja = 0.04. The dotted line is the xPT prediction for the 
pion and the sohd line is the prediction for the rho from IIJJ. Open symbols denote extrapolations to 

J ^ 0. 



the isovector 0~ and 1~ states, ie. the pion and rho, are consistent with the values 
found in [H |6]. Throughout the vacuum phase (i.e. ^ < /Xq) Mt, and Mp are more or 
less constant as expected for states with 5 = 0, although they do increase slightly but 
significantly for /ia > 0.2. For the pion this appears to be a j 7^ artifact vanishing in 
the j — >■ limit, but things are not so clear for the noisier rho. At and beyond onset at 
/ioO ~ 0.4 the pion and rho signals become much noiser as reflected in the error bars, 
but it is still possible to identify trends. The pion starts to become heavier at onset 
and appears to increase in mass monotonically with /i in the limit j — >^ 0. The j = 
xPT prediction M^, = M^{fi = 0)6'(/io — /i) + 2yu6'(/i — /io) [8] (dotted line) is followed 
in a qualitative sense. The increase of the pion mass post-onset is characteristic of a 
state formed from q and q with a symmetric combination of quantum numbers under 
the residual global symmetries (i.e. the Ps state in the notation of [8J) in a theory 
with Dyson index jSn = 1. The same trend is seen in QC2D simulations with staggered 
fermions in the adjoint representation [9J, but the opposite, namely a decrease in 
post-onset, is found in a theory with = 4 such as QC2D with fundamental staggered 
quarks 

By contrast post-onset the rho becomes significantly lighter, in agreement with the 
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result found in simulations on 4^ x 8 with significantly heavier quarks by the Hiroshima 
group [3j. This effect becomes stronger as j ^0. Reduction of Mp in a nuclear medium 
has been proposed to explain the low mass lepton pair enhancement observed in heavy 
ion collisions With the exception of the point at iia = 0.45 Mp seems to follow the 
predictions of an effective model description of spin-1 excitations in QC2D [H], shown 
as a solid line in Fig. [2l Our data suggest, however, that the transition is sharper than 
that predicted by the model. 

Results for isovector 0"*" and isoscalar 1+ states are omitted from Fig. [2] because the 
data are too noisy to fit, and for the isoscalar 1~ for the sake of clarity. The latter 
follows the rho almost exactly until /za > 0.5 at which point it becomes too noisy to 
measure. We do, however, include results for the isoscalar 0~. Since only the connected 
diagrams are currently being considered the difference between these two states and the 
pion and rho is just the sign of the anomalous term in (171191) . The fact that they are the 
same until onset shows that this term is negligible in this regime. Beyond onset, there 
is a small window in which the 0~ in the isoscalar channel is significantly lighter than 
the isovector, and is indeed roughly degenerate with the 1 = 00^ diquark of Fig. HI to 
be discussed below. 

The two remaining mesons shown in Fig. [21 the isovector 1+ and isoscalar 0+, show 
a similar behaviour, both starting off relatively heavy (and noisy) and then rapidly 
dropping as /i increases. By /ia = 0.3 they have reached a minimum and stay more 
or less constant as fi increases further. We shall argue below that the low mass of the 
0+ state is due to its overlap with the Goldstone boson in the superfiuid phase, when 
baryon number ceases to be a good quantum number; the low mass of the 1^ is more 
unexpected. 

To understand the physics of the diquark sector, it is helpful to begin with the Higgs 
and Goldstone states of Eq. (in|) with a varying diquark source strength j. Fig. 131 shows 
Higgs and Goldstone masses as functions of fi at ja = 0.04. The insets show how the 
two states scale with j at three selected values of /i. Below onset Higgs and Goldstone 
are degenerate, both scaling approximately linearly with j. Post onset the degeneracy 
is broken, and the relation Mcoid oc a/J predicted in xPT [8j appears to hold, the x^/^di 
being 0.67, 0.62 and 1.65 for the 3 linear Higgs extrapolations, and 1.08, 2.69 and 0.24 
for the 3 Goldstone extrapolations, which have the form Mcoid = ci^/j + ^ above onset. 

Extrapolations to j = using these scaling assumptions are also displayed on the 
main graph (triangular symbols). The two states remain degenerate until onset at which 
point the Goldstone becomes lighter than the Higgs, and appears to become massless 
as j —>■ 0. This is a clear manifestation of spontaneous breaking of U(1)b symmetry 
breaking for fi > fi^, implying a superfiuid ground state in which baryon number is no 
longer a good quantum number, and therefore meson and diquark states in principle 
indist inguishable . 

The diquark spectrum in the remaining spin-0 and spin-1 channels is shown in Fig. HI 

''The pre-onset behaviour M(j) = M{0){l + bf)i predicted by xPT [8] may be difficult to distingush 
from Unear behaviour in this regime. 
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Figure 3: Higgs and Goldstone masses as a function of fi. The two insets show results as j is varied at 
fixed iia = 0.3, 0.5, 0.7. Extrapolations to j = are displayed on the main graph with triangle symbols. 



It is striking that the signal-noise ratio is much higher for some diquarks than for the 
mesons, also seen in simulations with staggered fermions [1] . The two cleanest signals are 
for the isoscalar O"*" and the isovector l"*". The first observation is that as a consequence 
of the symmetry there is a relation between the meson and diquark spectra which 
holds for /i = j = if disconnected diagrams are neglected: 



MDiJ"") = Mm{J- 



(15) 



For < /X < /io, during which the physical ground state remains the vacuum, we thus 
predict Md{0^) = ± 2/i, Moil^) = Mp± 2/i, shown as dot-dashed lines in Fig. H 
Indeed both diquark particle-antiparticle pairs behave as expected up to /la ~ 0.3. 
Diquark masses are not shown beyond /xa = 0.25 as they become unfittable, as explained 
in Sec. [31 After this both O"*" and l"*" anti-diquark states flatten off and slowly decrease 
with /i. The other two isoscalar diquarks constructed from local operators, namely the 
0~ and 1~, are extremely heavy and hard to fit below onset, but above onset have a 
sufficiently good signal for us to deduce masses comparable with M^(yU = 0), Mp(/i = 0). 
Although the noise in the meson sector is admittedly large, the approximate degeneracy 
between meson and baryon sectors in the O"*" and 1^ channels seen in Figs. |2] and H] is 
consistent with the meson-baryon degeneracy in the superfiuid state discussed above. 
Meson-diquark degeneracy has also been observed in quenched studies at 7^ with 
staggered fermions [12] . 
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Figure 4: Diquark masses as a function of /i. Only states corresponding to operators with B < are 
displayed for /xa > 0.25. The dot-dashed hnes have intercepts at AI^^^p, and gradients ±2. 



The isoscalar 0~ diquark is a particularly interesting state in QC2D, since because 
of meson-baryon mixing in the superfluid phase it has the same quantum numbers as 
the 7]' meson [2j, and hence its mass acts as a probe of instanton effects and/or possible 
restoration of U(1)a symmetry in a baryonic medium. Unfortunately the current sim- 
ulations are not close enough to the chiral limit to settle this issue via observation of a 
7r(/x = 0)-?7'(/i) mass splitting. 

The kaon spectrum is shown in Fig. [5l since fig is assumed to be zero these states 
in effect carry a baryon charge conjugate to /i, and hence in general their correlators 
are not forwards-backwards symmetric. Some care must be taken in assigning physical 
quantum numbers to the states [5J. Consider say the 0~ excitation in a simple model in 
which the quarks are non-interacting but have mass m. Below onset the "light" quark 
spectrum is split into two branches, a hole branch with energy Eh = \/j^~-\-rn? — fi and a 
particle branch with Ep = \/p'^ + m? + /i, while the s-quark has energy Eg = a/p^ + 
independent of /i. We thus identify a zero momentum state with energy Eg + Eh = 

2m — /i, and a "ii'"'"" with Es + Ep = 2m + fi. Above onset, the particle branch becomes 
filled up to some momentum pf = ^/Jj^~-m^, so that the light quark hole state now 
costs energy Eh = fi — \/p^ + m?. The K~ state now has E = Es + Eh = n independent 
of p, implying a very small in-medium velocity. If the discussion is modified to include 
the effect of a non-zero diquark source, then the behaviour both below and above onset 
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1.25 




0.75 



Figure 5: Kaon spectrum as a function of /i. Only the K~ states are displayed for > 0.25. Dot- 
dashed lines correspond to Mk — Af^^p(/i = 0) + /i, while the dashed line corresponds to with 

p = and m = M-^^p{iJ. = 0)/2 



is described by the single solution 

E{p) = ^(yU — a/p^ _|_ m'iy -)- j2 _|_ ^p2 _|_ ^2 _|_ j2_ ^j^g-) 

Fig. [5] shows that this simple picture works qualitatively quite well in both K and K* 
channels, even once interactions are present, if we replace the vacuum mass 2m with the 
/i = mass of the corresponding meson from Fig. [2J The observed Mx always exceeds 
the model prediction, but appears to match both at = and /i — >• cxo and retains a 
minimum for /i ^ /io. There is thus a significant range /Xo < /i ^ Mk in which M^- 
lies below its vacuum value, and similarly for Mk*-- This offers theoretical support 
for the idea that in-medium modification of kaon masses in nuclear matter leads to 
deeply-bound /^"-nuclear states [TH], for which there is some experimental support from 
scattering kaons on light nuclei [H]. Of course, simulations with realistic quark masses 
are needed before this observation can be taken seriously, not to mention the extension 
from two colours to three. Nonetheless, we can identify a minimal requirement for the 
existence of kaonic-nuclear bound states, namely Mk > fJ-o- This inequality is satisfied 
in QCD. 
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Quantum Numbers 




Meson 




Diquark 








/i < /io 


/i > /Xo 


/i < /io 


/i > /io 


/ = 


0+ 


noisy 


good 


good 


good 


/ = 


0- 


good 


noisy 


noisy 


good 


/ = 




good 


noisy 


noisy 


good 


/ = 


1+ 


no signal 


no signal 


no state 




/ = 1 


0+ 


no signal 


no signal 


no state 




/ = 1 


0- 


good 


noisy 


no state 




/ = 1 


1- 


good 


noisy 


no state 




/ = 1 


1+ 


noisy 


good 


good 


good 



Table 2: Informal summary of the success of our hadron fits. 



5 Discussion 

It is helpful to summarise our findings in the form of a table in which the quality of 
our fits both below and above onset is given. The main achievements of this study have 
been observations of: 

• The reversal of the pion and rho levels on crossing from vacuum into a baryonic 
medium. In the vacuum /i < /io M^- p is approximately constant, probably because 
there is no diquark state with the same quantum numbers with which to mix. 

• The breaking of the degeneracy between Higgs and Goldstone diquark states for 
/i > /io, the Goldstone mass scaling as \f] in accordance with general theoretical 
properties of spontaneous symmetry breaking by condensation of fermion pairs. 

• Further evidence for meson-baryon mixing in the degeneracy of / = 0^ and / = 1 
1"*" states for /i > /io. Post onset the l"*" appears to be the next lightest state after 
the Goldstone and Higgs. The fact that the mesons with these quantum numbers 
appear not to have constant mass even pre-onset (see Fig. [2]) can also be ascribed 
to meson-baryon mixing, since for j 7^ there is a non-zero amplitude for ■?/''?/' to 
project onto a baryon. 

• Evidence of the possibility of bound kaons in nuclear matter. 

Our study has uncovered little about the effects of a second deconfining phase tran- 
sition suspected to occur at /i^a ^ 0.65 on this system [1]. The only possible discernable 
trend is a levelling off of the already massive pion state for /xa ^ 0.5 seen in Fig. [2J Sta- 
tistical noise increasing with /i, however, makes this observation provisional at best. In 
a deconfined phase we might expect mesons and baryons to be formed from particle-hole 
and particle-particle pairs in the neighbourhood of a Fermi surface, and it is possible 
that the local operators used in this study have a poor projection onto the true quasipar- 
ticle excitations. We hope that studies of meson and diquark wavefunctions currently 
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in progress, as well as extended operators for the 7 = 00+ and 7=11+ mesons and 
enhanced statistics for the disconnected diagrams, will clarify the situation. 
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